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Abstract 

It has been argued recently that the inclusion of the chiral logarithms in ex¬ 
trapolation of the lattice data can shift the value of the hadronic parameter £ = 

f Bs \J'Bb s / fB d \]'RB d , from 1.16(6) to 1.32(10) and even higher. If true, that would 
considerably change the theoretical estimate for the ratio of oscillation frequencies 
in the B®- and -systems, and would affect the standard CKM unitarity triangle 
analysis. In this letter we show that f Bs / f Bd ~ fx/fm and thus the uncertainty due 
to the missing chiral logs is smaller than previously thought. By combining the NLO 
chiral expansion with the static heavy quark limit we obtain £ = 1.22(8). 


1 Uncertainty in comes from fsjfB d 

The K° — K° and B°—B° mixing amplitudes are the main ingredients in the standard unitar¬ 
ity triangle analysis [I] . They overconstrain the apex of the CKM triangle in the (p, r))-plane. 
Confronting such an overconstrained triangle with the directly measured sin(2/3) [5] allows 
for the consistency check from which we hope to either see the effects of the non-Standard 
Model physics, or simply to confirm the validity of Standard Model at the accessible accu¬ 
racy J3j- One of the essential constraints is provided by the ratio of oscillation frequencies 
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where rn Bs(d) is the mass of meson, while Vt s and Vtd are the CKM matrix elements. 
The mass difference A m,B d , which measures the oscillation frequency in the B d — B® system, 
has been determined accurately in the experiments, Am Bd = 0.503 ±0.006 ps _1 |4]. On the 
other hand, Amg s is currently only bounded from below, A m Bs > 14.4 ps -1 [4J, but will 
hopefully be measured soon at Tevatron ji|. On the r.li.s. of eq. (HD, one has the hadronic 
parameter 
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with f Bq and B Bq being the decay constant and the so called “bag”-parameter, respectively. 
The value of the parameter £ is expected to be determined from the QCD simulations on 
the lattice with a good accuracy since many uncertainties should cancel in the ratio. 

Recently, however, it has been argued that the chiral logarithms might be a source of 
a large, yet unaccounted for, uncertainty in £ 0|. The reason for this lies in the fact that 
the masses of light quarks that are directly accessible from the lattice studies are those 
around the strange quark, m s /2 < m q < 3m s /2. In this range f Bq and B Bq , computed on 
the lattice, exhibit a linear behaviour under the variation of m q . This linear dependence 
is then extrapolated in m q —» m^, to reach f Bd and B Bd . The claim of ref. 0 is that the 
inclusion of the chiral logarithms in the mentioned extrapolation produces a shift from the 
“standard” value f B Jf Bd = 1.16(5), to 1.32(8). In ref. J7J the systematic error due to the 
shift is estimated to be even larger, ±0.24. In the following we revise the problem and argue 
that the shift is indeed present but it is actually smaller than previously claimed. We will 
show that f B JfB d ~ fx/fn- 

Finally, from refs. mm it is known that, w.r.t. the result of the linear extrapolation 
in the light quark mass, the chiral logs produce a tiny (insignificant) shift in the ratio of 
“bag”-parameters, B b JB B d . This latter ratio is known to be very consistent with unity ||, 
and therefore is not the source of any additional uncertainty in the determination of £. 


2 Ratios fn s /fn d and fi±/fn from the ChPT 


The SU(3) light flavour breaking effects in the ratios of the decay constants have been 
calculated at the 1-loop level of the chiral perturbation theory (ChPT), resulting in the 
following expressions mmxn 
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with the convention / = 130 MeV, and the standard notation $ Bq = f Bq y/m Bq . The 
counter-term L 5 is defined in im, while K stands for a sum of counter-terms discussed in 


2 










im Note that these counter-terms provide the corrections to linear in quark mass 

and that they eliminate the renormalization scale dependence arising from the chiral loop 
integral 


I\{m 2 ) = m 2 log(m 2 //i 2 ) . (5) 

The renormalization scheme of ref. cm has been used. The only new parameter accompany¬ 
ing the chiral logarithm in eq. © is g, the coupling of the lowest lying spin doublet of heavy 
mesons to a soft pion. Its experimental value has been recently established in the charm sec¬ 
tor, namely g c = 0.59(8) [12] • While a recent (quenched) lattice study suggests that the g- 
coupling does not change when increasing the heavy quark mass (g b ~ g c ) ^, the light cone 
QCD sum rule calculation predicts the suppression by about 25% (i.e. gb/g c ~ 0.75) T3| . 
We will take the average of the two predictions and add the difference in the systematic 
uncertainty, i.e. g = g b = 0.52(7)(7). It is worth mentioning that the use of the tree 
level values for g and / is consistent with the use of the ratios © and ©^ derived at 
next-to-leading order in the chiral expansion. 

Eqs. © and © have been obtained in refs. nu and cm. respectively. Notice that in the 
former case we neglect the OQ/mQ-corrections, which are anyway expected to cancel in 
the ratio © to a large extent HU. Throughout this letter we assume the isospin symmetry 
rrid = m u = rn q , and use the Gell-Mann-Oakes-Renner and Gell-Mann-0kubo formulas to 
write 

2 2 r + 1 2 'r-j-2 . . 

m n = 2 B 0 m s r , m~ K = 2 B 0 m s —^— , m n = 2 B 0 m s —^— , (6) 

where r = m q /m s: and B 0 is the common factor related to the quark condensate, B 0 = 
—2 (qq)/f 2 - When varying the ratio r, we will always keep the strange quark mass fixed 
to its physical value. We can then use the physical kaon and pion masses to get 2 B 0 m s = 
2(m T £' ys ) 2 — ( m P hys ) 2 = 0.468 GeV 2 . 

It is important to note that even though the chiral logarithms in eqs. © and © enter 
with different coefficients, they are numerically very similar. To illustrate this point we set 
K = L 5 = 0 and plot in fig. ID only the chiral log contributions to Rf B and Rf n . Indeed, we 
see that for the central value of the parameter g, the chiral log contributions to Rj b and 
Rf n nearly coincide, whereas the relative difference at the physical point r ~ 0.04 is below 
5% level when g is varied in the range g G (0.35,0.7). Note that the near coincidence of 
chiral logarithmic terms in Rf B and Rf n is not spoiled by the change of the renormalization 
scale /r. 

From fig. HI we also observe the well known fact that the chiral logs are not large enough 
to reproduce the experimental value (Jr/ fir)ex P . = 1.22(1) HU- The mismatch is patched 
by fixing the coupling = 1 GeV) to (8.7 ±2.5) x 10~ 4 from eq. ©. The remaining 
entry in eq. ©, the value of K, will be discussed later on, at the end of section H 

3 Chiral extrapolation of fB s /fB q 

Contrary to the ChPT, which is valid for very small quark masses, the lattice QCD provides 
predictions for not so light quarks corresponding to 1/2 < r < 3/2. In other words, fs a 
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Figure 1: The chiral logarithmic contributions to Rf„ and Rf B for three different values of g spanning the 
range g = 0.52(7)(7), discussed in the text. Plotted are eqs. © and © with K = L 5 = 0 and g = 1 GeV. 
Variation of r = m q /m s is made by keeping m s fixed to its physical value. 


(r = r s = 1) is accessed directly from the lattice studies, while for fs d an extrapolation in 
r —> rd is needed (r^ ~ 0.04 [03). That extrapolation is typically made linearly as 

i2jf = l + a(l-r) , (7) 

with a = 0.16(5) jj, I?;. The form © does not include the chiral logarithmic terms shown 
in eq. ©, which would increase the result of extrapolation, namely fs a / fB d - It is, however, 
not clear at what value of r the leading chiral logs need to be included, i.e. we do not know 
how light the “pions” should be so that the leading chiral logs remain uncompensated by 
the higher order terms in the chiral expansion. 

In the approach adopted in ref. [6] the NLO prediction of ChPT © is assumed to be 
valid up to r = 1, while the counter-term K is fixed from the requirement 
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at some intermediate point r^r- We slightly modify that approach by actually separating 
2 regions: (i) above r M , where the extrapolation of the linearly fitted lattice data © is a 
good approximation; (ii) below tm, where we use the ChPT formula © with K fixed from 
the condition ©, and go down to the physical point The effect is shown in fig. ©where 
we plot two curves, corresponding to tm = 0.25, and 0.75 1 . From fig. ©we see that farther 

1 For an easier orientation we convert r« by using rn' 2 M = 2Bom s rM (2Bom s = 0.468 GeV 2 ) to obtain 

rrijvf « 330 and 580 MeV, respectively. 
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Figure 2: Chiral extrapolation of the ratio Rj b = fsj $B q , in which the chiral logarithms © are included 
from tm = 0.25 and 0.75 with K fixed from © as in ref. [5]. Note that the variation of r = m q /m s is 
made by keeping m s fixed to its physical value. The illustration is provided for g = 0.52. We also show 
Rf B = /bJ fB q with K fixed as I\ = L 5 (dotted line). 


away from lies the point rjf, the larger is the shift in Rj b - By using g = 0.52, we obtain 2 

Rf B = l-22(6) rM=0 .25, 1.30(7) rM= o. 75 , (9) 

to be compared with the result of the linear extrapolation Rj b = 1.16(5). To estimate this 
shift, the authors of ref. [0] take = 0.75 ± 0.25. Although the relevance of the leading 
chiral logs for the matching scale at rriM — (580 ± 100) MeV may be doubted, at this point 
there is no legitimate argument to exclude the possibility that the central value for Rf B 
gets as high as 1.3. 


4 Double ratio (f B JfB d )/(,fi</ A) 

To get around the above ambiguity we construct the double ratio [20] 
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2 Had we used the approach similar to the one in ref. 03 with both the size and the derivative of Rf B 
matched to the lattice data at Tm , the corresponding shifts would be smaller than the ones presented in ©. 
Note that this procedure amounts to introducing an additional constant term K' in © (arising from NNLO 
contributions) with both K and K' then fixed from lattice. 
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As shown at the end of section [21 the chiral log contributions in R = Rf B /Rf n almost 
completely cancel so that the double ratio R , to a very good approximation, depends 
on r only linearly. One can then envisage a future lattice analysis in which the chiral 
extrapolation of R is made linearly and then the result for fs 3 / fB d is deduced by using the 
experimental value [Ik/ In) exp. = 1.22(1). 

In this letter, we will use the formula era to estimate the value of /s s //s d . To do so 
we have to fix the couplings L 5 and K. We already mentioned that L 5 = (8.7 ± 2.5) x ICO 4 
(at fi = 1 GeV). However, the coupling K is unknown. To get its value, but only relative 
to L5, we may rely on the experience with the QCD simulations on the lattice in which the 
heavy-light (B q ) and light-light (P qq ) pseudoscalar meson decay constants are fitted well to 
the linear forms 


fB q = a 0 + a\rri 2 Pqq , and f Pqql = b 0 + . 

By using eq. these forms can be rewritten as 


R fg ~ 1 + — 2B 0 m s (l - r ) , 
Po 
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The values of parameters a and (3 can be obtained from the linearly fitted lattice data for 
Rf B and Rf n , which at r = 0 simply read, R.f B = 1 + a and Rf n = 1 + (3. From the recent 
compilation of unquenched calculations (with rif = 2 ) of the heavy-light decay constants, 
which are obtained through the linear extrapolation, we take fnjfn d = 1.16(5) 0 T7j . 
The MILC collaboration has recently verified that the ratio /e s //s d remains stable when 
the number of dynamical flavors rif and/or the values of their masses are changed (HU) . 
We therefore take a = 0.16(5). From an extensive study of the light hadrons with two 
flavors of dynamical quarks by the JLQCD and the CP-PACS collaborations [21], one has 
f k/ fn = 1.16(2), i.e. f3 = 0.16(2). 3 I11 other words, a//3 ~ 1. That conclusion is also in 
agreement with the observation made by the MILC collaboration that the slope bi is about 
30% larger than the slope Gq [22] 4 - Moreover, by fitting, linearly in r, the unquenched 
data for both the light-light and the heavy-light decay constants obtained by the JLQCD 
collaboration and presented in ref. one arrives to the same conclusion, namely a/(3 ~ 1. 

This information can be used to estimate the value of K — L5 by matching at some 
intermediate point the chiral expression for R given in eq. m to the linear fits of 
lattice data given in eq. m 
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In further discussion we will take to be in the range Tm £ (0.5,1), where lattice simu¬ 
lations are performed. From requirement m it follows trivially 


of p hys\2 1 

_ - L f J ~(K ~ In) = (<*-/»- M’-U,,, ■ (14) 

3 This result corresponds to the average of the results obtained by using K- and 4>~ inputs in the recent 
paper by the JLQCD Collaboration in which the results of the improved unquenched study were reported. 
4 We thank Claude Bernard for this information. 
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where 5R loop is the chiral log correction to the R ratio (i.e. the first correction term in 
eq. (□m while a and /3 are the slopes of the fits to the lattice data points as defined in 
eq. EH- Finally, the expression for the double ratio eh is 

R = 1 + 5R loop +((a-P) --— 5R loop 

V 1 -r M 

As we already mentioned, the available lattice data suggest a ~ (3. Furthermore, the size 
of chiral corrections to the matching condition, d-R loop | r _, r /(I — r^f), is below 1 % for 
r M E (0.5,1) and g E (0.35, 0.7). As such it is much smaller then the present errors on the 
value of a — (3 and can be safely neglected in the following. 

As a first estimate we can thus set the last term in eq. m to zero. At the physical 
point r = rd = 0.04, and by using g = 0.52(7)(7), we get 

R = 1.011 ± 0.024 . (16) 

where the errors reflect the variation of the (^-parameter only. From this result, we conclude 

that f Bs /fB d ~ fx/fn- 


I r=r M 


(1 - r) 


(15) 


5 Our estimate of fnjfn d and £ 


Finally, to make the numerical estimate we take 

/ = 0.13 GeV, 0 = O.52(7)(7), K/L 5 = 1.0 ± 0.6 , L 5 = (8.7 ± 2.5) x 10 “ 4 , (17) 

where the first error in g is considered to be gaussian (experimental) and all the others to 
be flat. The error in (fx/fPexp. = 1.22(1) is also gaussian. Notice in particular that the 
value of K/L 5 ratio follows from the present value of a — f3 = 0.00 ± 0.06, discussed in the 
previous section, through the use of eq. EH- That uncertainty may easily be reduced by 
the direct lattice calculation of the ratio R. 

With these input values and by using eq. EH, we produced the histogram of 10 6 Monte 
Carlo events for f B Jf Bd = R\f^ B J^ Bs x ( Ik/ fn)ex P ., which is shown in fig. El We finally 
obtain 


f B Jf Bd = 1.22 ±0.08, 


(18) 


which is larger than the commonly quoted values, yet it is considerably smaller than the 
results obtained in refs. 0 Gl¬ 
ib make the discussion complete and get an estimate of the parameter £, we will also 
use the result of ref. 
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where K denotes the sum of corresponding counter-terms. Combining eq. EH with eq. EH 
gives 
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Figure 3: Probability distribution function for the values of the SU(3) breaking ratio fB s /fB d , as obtained 
by using the input values discussed in the text (see eq. 1171) 'l. 


As mentioned at the beginning of this letter, the results of several lattice studies suggest 
that K can be safely set to zero j9|. With the same input values, given in eq. , we 
compute lZ£^/mB d /rnB s x (fx/fir)ex P . from eq. (J20J), and obtain 

f = 1.22 ±0.08. (21) 

ft is important to stress that our results for JeJIr,, and £ are obtained by neglecting 
the O(l/m^ >0 ) corrections that survive the cancellation in the ratio (J3J). One may worry 
about the possibility that those corrections may further increase our results. The experience 
with lattice QCD and with the QCD sum rule (QSR) calculations suggests that the SU(3) 
breaking ratio of the decay constants is not decreasing when increasing the heavy meson 
mass. Indeed, from ref. H7| we read the lattice QCD estimates, / b s /Ib = 1.16(5), fnjfo = 
1.12(4), while from the last paper of ref. jOJ we see that the QSR calculation gives, /s 3 //b = 
1.16(5), fDs/f d = 1.15(4). Therefore the uncancelled corrections arc likely to shift 

our result m towards the smaller values rather than the larger ones. 

Notice that a conclusion similar to ours’ has been also reached in the recent reanalysis 
of the lattice data produced by the MILC collaboration [22] ■ The authors show that the 
consistent inclusion of the chiral logarithms in all the chiral extrapolations involved in their 
data analysis, produces a shift of their central result, fsjfB d = 1-16(1)(2)(2), by only 
+0.04, which is much smaller than +0.24, as claimed in ref. [7], or +0.16, as argued in 
ref. |L Moreover, by actually shifting their result, one gets fsJfB d = 1.20(1)(2)(2), which 
is completely consistent with our main claim, namely fB s /fB d ~ fi</ fir- 

Finally we mention, that instead of the double ratio ( fB s /fB d )/(fK/fn ) discussed in the 
present letter, the double ratio of the heavy meson decay constants Ri = ( fB s /fB d )/(fD s /fD d ) 
could be used for the determination of (/b s / fB d ) and correspondingly of £ as suggested by 



































Grinstein in [21] . The ratio of charm meson decay constants (/d s //dJ is expected to be 
measured by CLEO-c, while preliminary Nf = 2 dynamical lattice calculations of R\ are 
already available f? . 

6 Conclusions 

In this letter we have discussed the effect of the inclusion of the chiral logarithms on the 
lattice determination of the ratio fn s /fn A ■ We have shown that in the double ratio 

j_) fB s / f B d 

the chiral logarithms almost completely cancel because of the fact that the numerical value 
of the coupling g is close to ~ 0.5. To a very good approximation the double ratio R then 
depends linearly on the light quark mass all the way down to the chiral limit. 

We have also made a numerical estimate of the hadronic parameter £, by using the 
expressions derived in the ChPT and by fixing the unknown couplings from the available 
lattice results. Namely, we obtain 


f = 1.22 ±0.08. 
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